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1.  INTRODUCTION 

This  paper  develops  models  of  the  transient  behavior  of  Markovian 
repairable  item  provisioning  systems.  A  multi-echelon  structure  of  re¬ 
pair  and  resupply  is  of  concern.  It  is  desired  to  develop  analytical 
solution  techniques  for  exact  models  of  systems  with  finite  numbers  of 
items  and  finite  repair  capacities  (i.e.,  a  finite  number  of  repair 
channels  at  each  repair  facility) .  This  is  in  contrast  to  Dyna-METRIC 
[see  Hlllestad  (1981)  and  Hillestad  and  Carrillo  (1980)],  which  is  an 
approximate  model  assuming  an  infinite  calling  population  of  items  and 
ample  repair  capacities. 

Most  multi-echelon  repairable  item  work  has  concentrated  on 
steady  state  solutions  and  revolves  around  the  METRIC  type  of  model, 
which  assumes  an  infinite  population  of  items  which  can  fail  and  an 
infinite  repair  capacity,  so  that  no  queue  ever  forms  at  the  repair 


« 
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facility  [see  Sherbrooke  (1968)  and  Muckstadt  (1973),  for  example]. 
Recently,  the  METRIC  type  of  model  has  been  extended  to  provide  tran¬ 
sient  solutions  for  a  time  varying  environment  (Hillestad  and  Hillestad 
and  Carrillo,  op.  oit.) ,  and  is  called  Dyna-METRIC. 

Exact  models  for  finite  item  populations  and  repair  capacities 
have  been  studied  by  Gross,  Miller,  and  Soland  (1981),  but  for  steady 
state  solutions  only.  This  paper  presents  results  on  transient  solu¬ 
tions  for  such  problems. 

For  convenience  in  describing  the  various  multi-echelon  systems 
studied,  the  following  classification  scheme  is  adopted.  A  system  is 
described  by  seven  symbols,  in  the  format  (n^.n^.n^/n^.n^/ng.n^) . 

Table  1  provides  the  definitions  of  n^  through  n^. 

As  an  example  of  the  classification  scheme,  consider  a  military 
supply  system  with  five  bases,  base  and  depot  repair,  base  and  depot 

TABLE  I 

SYSTEM  CLASSIFICATION  SCHEME 

Symbol  Definition 

n^  Number  of  bases  (1,2,...) 

n^  Number  of  levels  of  repair  (1,2,...) 

n^  Number  of  levels  of  supply  (1,2,...) 

n^  Number  of  item  types  (1,2,...) 

n^  Number  of  levels  of  indenture  (0,1,...) 

n^  Size  of  item  population  (f  or  ®) 

tiy  Capacity  of  repair  facilities  (f  or  «) 
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•paras,  two  types  of  Items  with  no  Indenture,  finite  calling  population 
slse  (finite  numbers  of  Items),  and  finite  repair  capacities.  This 
system  would  be  represented  as  (5,2,2/2,0/f ,f) .  The  same  system  with 
Infinite  calling  population  and  infinite  repair  capacity  (Dyna-METRIC 
aseiaptlons)  would  be  denoted  (5, 2, 2/2,0/*,*) .  The  most  general  system 
descriptor  would  be  (b,r,s/k,J/f ,f) ;  that  Is,  a  system  with  b  bases, 
r  levels  of  repair,  s  levels  of  supply,  k  item  types,  j  levels  of 
Indenture,  and  finite  calling  populations  and  repair  capacities.  To 
describe  fully  any  system,  it  would  still  be  necessary  actually  to  draw 
the  network  structure;  nevertheless,  this  classification  scheme  will  be 
most  useful  in  delineating  the  systems  under  discussion. 

Our  goal  is  to  present  exact  models  of  Markovian  repairable  item 
provisioning  systems  which  can  be  solved  numerically  for  exact  transient 
state  probabilities  and  other  performance  measures.  We  are  also  con¬ 
cerned  with  the  practical  problem  of  algorithmic  implementation  on  a 
computer.  It  is  possible  to  solve  nontrivial,  fairly  complex,  systems 
using  the  algorithmic  approach  of  "randomization."  We  have  developed 
highly  efficient  algorithms  for  two  cases  [ (1,1,1/1,0/f ,f)  and  (2,2,2/ 
l,0/f,f)],  and  discuss  generalizations  of  these  situations. 

The  approach  we  use  is  a  two  stage  procedure.  The  first  stage 
is  a  method  of  modeling  called  SERT  and  the  second  is  an  algorithm  based 
on  the  randomization  technique  [the  reader  is  referred  to  Gross  and 
Miller  (1982)  for  a  full  methodological  discussion  of  these] .  A  good 
discussion  of  the  advantages  of  the  randomization  procedure  as  a  numeri¬ 
cal  analysis  tool  is  found  in  Grassmann  (1977). 
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2.  SERT  MODELING  APPROACH 

Since  we  assume  a  finite  Markovian  system  (all  failure  times  and 
repair  times  are  independent  exponential  random  variables) ,  the  system 
can  be  fully  described  by  its  infinitesimal  generator  matrix 

4911. 

it-*>  u 

where  P(At)  *  (p^(At)};  p^  (At)  =  Pr{system  is  in  state  j  at  time 
t  +  At  |  in  state  i  at  time  t}  .  The  matrix  Q  is  finite,  say  N  *  N  , 
since  we  are  assuming  a  finite  number  of  items  in  the  system.  For 
complex  systems,  N  can  be  quite  large  (this  will  be  seen  later),  so 
it  is  necessary  to  have  algorithms  as  efficient  as  possible,  both  in 
running  times  and  storage  demands  made  on  the  computer. 

It  is  not  necessary  to  store  this  N  *  N  Q  matrix  if  we  use 
the  SERF  modeling  approach.  Briefly  [for  more  detail  see  Gross  and 
Miller,  op.  ait.],  we  must  describe  the  State  space,  the  types  of  Events, 
the  transition  Rates  (the  off-diagonal  nonzero  elements  of  the  Q  matrix), 
and  the  Target  states,  that  is,  the  state  to  which  the  system  goes  when  a 
given  type  of  event  occurs. 

Given  a  state  space  S  of  size  N  with  s  denoting  a  given 

state  of  the  system  (s  e  S) ,  and  an  event  space  E  of  size  E  with 

e^  denoting  an  event  of  type  J  (e^  e  E) ,  it  is  necessary  to  consider 

a  rate  vector  and  a  target  state  vector  for  each  e^  ,  which  we  denote 

by  r^  and  t^  ,  respectively.  The  dimensions  of  r^  and  t^  are 

1  x  N  .  Since  in  these  models  E  «  N  ,  it  is  much  more  efficient  with 

respect  to  computer  storage  requirements  to  work  only  with  the  r^  and 

t^  vectors,  rather  than  the  N  x  N  Q  matrix,  as  the  former  requires 

2 

consideration  of  only  2*E*N  elements  rather  than  N  . 
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The  SERT  modeling  procedure  will  be  illustrated  on  the  (1,1, 1/1,0/ 
f,f)  and  (2,2,2/1,0/f ,f)  systems  to  follow. 


3.  RANDOMIZATION  COMPUTATIONAL  TECHNIQUE 

Transient  solutions  to  Markovian  queues  require  solving  the  set 


of  differential  equations 


ir'(t)  -  ff(t)Q 


where  ir(t)  is  the  vector  of  transient  probabilities,  i.e.,  ir(t)  ■ 

{7T .  (t) }  ,  where  n.  (t)  =  Pr(system  is  in  state  j  at  time  t}  and  ir'(t) 
j  3  ”* 

is  the  vector  of  derivatives  with  respect  to  time,  i.e.,  ir’(t)  *  (ir'(t)}  . 

-  j 

Many  numerical  techniques  can  be  used  to  solve  these  linear,  first 
order  differential  equations;  for  example,  numerical  integration  [see 
Maron  (1982) ] .  We  choose  instead  a  technique  referred  to  as  randomiza¬ 
tion  [see  Gross  and  Miller,  op.  ait.,  or  Grassmann  (1977)],  which  is 
Ideally  suited  to  the  SERT  modeling  of  Markovian  systems. 

Randomization  is  based  upon  the  ability  to  transform  the  continu¬ 
ous  parameter  Markov  process  analysis  to  an  analysis  of  a  discrete  pa¬ 
rameter  Markov  chain  (MC)  whose  transition  times  are  generated  by  a 
Poisson  process. 

Consider  the  continuous  parameter  process  we  are  modeling.  It 
remains  in  a  given  state  s  until  one  of  many  (E  possible)  events  occurs, 
which  then  changes  its  state  to,  say,  s'  (if  it  is  an  event  of  type 
ej  ,  s'  *  t^).  Since  all  times  are  exponential,  the  time  to  the  next 
event  is  the  minimum  of  exponentials  which  is  also  exponential.  Hence, 
if  we  look  at  the  process  only  at  transition  times,  it  is  a  discrete 
psrameter  Markov  chain,  whose  holding  (transition)  times  are  exponentially 
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Since  the  mean  of  the  holding  time  in  state  s  may  depend  on  s  , 
we  do  not  quite  have  a  Poisson  process  generating  these  transitions.  To 
get  around  this,  we  denote  the  mean  of  the  minimum  holding  time  over  all 
states  (this  always  exists  since  S  is  finite)  as  1/A  and  consider  a 
Poisson  process  with  rate  A  as  a  generator  of  the  transitions.  Be¬ 
cause  this  generates  transitions  at  a  greater  rate  than  desired  (since  A 
is  the  maximum  of  the  state-dependent  votes) ,  we  must  thin  the  process 
to  model  the  actual  state-dependent  holding  times  by  adjusting  the 
discrete  parameter  HC  transition  probability  matrix.  Without  the  ad¬ 
justment,  the  embedded  discrete  parameter  MC  has  transition  probability 


matrix  P  ■  {p  . }  ,  p, .  ■  q. ./  £  q..  .  The  parameter  A  is 
max  £  q^j  .  Adjusting  the  P  matrix  for  the  thinning  operation  gives 


i  j*i 


the  transition  probability  matrix  P  *  {p^}  *  where 


l  q,.,  Q 

liJ _ _  ,  hi 

'ij  7  A  A  * 

jl 


_ * 


Thus,  the  computing  formula  for  ir(t)  ■  {it  (t) }  is 

—  s 


<*>  n 

*3<t)  *  l  l  Mn> 

n-0  i-0 


e~At(At)n 

n! 


(i) 


where  <j>  (n)  is  the  probability  that  the  system  is  in  state  s  after 
n  transitions  generated  by  the  Poisson  (A)  process,  and  e  At(At)n/n! 
is  the  probability  that  there  are  n  transitions  of  the  Poisson  (A) 
process  in  time  (0,t)  .  The  $8(n)  can  be  determined  using  ?  and 
it^(O)  in  the  standard  MC  way;  that  is, 

♦  <n)  -  ir(0)(f)n 
*  <Kn  -  1)P  , 
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where  tt(0)  is  the  starting  state  probability  vector.  The  infinite  sum 
in  (1)  must  be  truncated  at  some  appropriate  point  (see  Gross  and  Miller, 
op.  cit.)  which  can  be  set  to  guarantee  a  bound  on  the  error. 

Using  the  target  and  rate  vectors  from  the  SECT  procedure,  we  can 
compute  the  J>(n)  vectors  recursively  in  a  manner  more  efficient  than 
using  (2)  by  the  following  algorithm: 

(i)  <>«»  -  ir(0) 

(ii)  <j>(n  +  1)  is  computed  from  4>(n)  as  follows: 

(a)  4>s (n  +  1)  »  ♦gin)  * 

then 


1  - 


E 

l 


(3) 


(b) 


for  j 

Vn) 


»  1 , 2 , . . . ,  E  ,  and  s  £  S  , 

•  (r^/A)  to  <j>  (n  +  1)  . 
s  tJ 


add 


s 

What  this  algorithm  does  is  to  operate  simultaneously  on  all  com¬ 
ponents  of  <}>(n  +  1)  using  the  components  of  <p(n )  and  the  transition 

probabilities  from  P  .  Note  that  (r^/A)  is  p  ,  ;  that  is,  the 

S  s,t^ 

probability  of  going  from  state  s  to  state  t^  s  (the  target  state 

that  event  j  causes  the  system  to  switch  to  when  it  is  in  state  s  ) 

given  an  occurrence  of  the  Poisson  (A)  process.  The  term 
E  1 

(1  -  .  rJ  /  A)  is  the  probability  of  ignoring  an  occurrence  of  the 

Poisson  (A)  process,  which  is  called  thinning  and  which  leaves  the  sys¬ 
tem  in  state  j  (a  null  event). 

This  algorithm  is  very  efficient  for  sparse  Q  matrices,  which  is 
the  situation  we  generally  have  for  repairable  item  systems,  since  it 
avoids  the  zero  multiplication  that  would  result  from  using  the  matrix 
multiplication  of  (2) . 
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4.  SIMPLE  NETWORK  MODEL:  (1,1,1/1,0/f ,f)  SYSTEM 

The  simplest  case  of  a  repairable  item  system  is  the  simple  ma¬ 
chine  repair  model  shown  in  Figure  1.  The  situation  modeled  has  a  popu¬ 
lation  consisting  of  M  items  desired  to  be  operational  at  all  times  and 
Y  spares  to  support  the  system.  There  are  c  repair  channels,  so  that 
a  maximum  of  c  items  can  be  undergoing  repair  simultaneously.  If  more 
than  c  items  require  repair,  a  queue  forms  at  the  repair  facility. 
Failure  and  repair  times  are  exponentially  distributed  random  variables 
with  the  mean  time  to  failure  of  any  item  denoted  by  1/X  and  the  mean 
time  to  repair  an  item  denoted  by  1/y  ;  i.e.,  X  and  y  are  the  failure 
rate  and  repair  rate,  respectively.  The  total  number  of  items  in  the  sys¬ 
tem  is  M  +  Y  =  N  and  if  s  items  are  in  repair,  then  N  -  s  items 
are  at  the  operating  node,  so  that  when  s  >  Y  ,  the  population  is  oper¬ 
ating  at  less  than  the  desired  strength  of  M  . 

For  modeling  purposes,  the  operating  node  of  the  two-node  network 
pictured  in  Figure  1  can  be  considered  to  be  an  M  channel  queue,  so 


z 

hi 


Figure  1. — (1,1,1/1,0/f ,f)  system. 
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that  if  there  are  more  than  M  units  at  this  node  (s  <  Y) ,  the  queue 
represents  on-hand  spares.  If  no  queue  is  present  and  some  "servers" 
are  idle  at  the  operating  node  (s  >  Y) ,  then  the  population  is  operating 
at  degraded  strength  (as  mentioned  above) . 

Utilizing  the  SERT  methodology,  the  state  of  the  system  s  can 
be  considered  to  be  the  number  of  items  in  or  awaiting  repair  (in  re¬ 
supply).  The  size  of  S  is  N  ■  M  +  Y  +  1  ,  since  s  ■  0,1,2 . M  +  Y 

There  are  only  two  types  of  events  (E  ■  2) .  The  first  type  of 

event  is  a  failure.  The  rate  vector  for  failures  is  rf  -  (r^,  r^,  .... 

f  f. 

rg,  ....  rN)  ,  where 


{MX 

(M  +  Y  - 


and  the  target  state  vector  is  t 

f 


!s  +  1 
N 


0  <  s  <  Y 
OX  Y  <  s  N 

-  tJ  J  -f 

\£q*  ...,  tg,  ...,  t^)  ,  where 

s  »  0,1, ... ,N  -  1 
s  ■  N 


since  a  failure  increases  the  number  in  resupply  by  one,  except  when  N 
are  in  resupply.  In  this  case  a  failure  cannot  occur,  but  we  show  t^ 
as  N  ;  that  is,  a  "mythical"  failure  does  not  change  the  state. 

The  second  type  of  event  is  a  service  completion.  The  sth  com¬ 
ponent  of  the  rate  vector  rr  is 

[sy 


0  <  s  <  c 

cy  c  <  s  <  N 

and  the  corresponding  sth  component  of  the  tr  vector  is 


r  (s-1  1  <  s  <  N 

8  io  s  -  0  . 
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We  can  now  apply  the  SERT  multiplication  algorithm  (3) ,  which  reduces  to 

(i)  4>(0)  -  TT (0) 

• ■ ^ 

(ii)  (a)  For  s  »  0,1,..., N  ,  set 


then 

(bL) 


$_(n  +  1)  ■  $  (n) 
s  s 


For  s  ■  0,1,..., N  ,  add 

f 

r 

♦8<»)  *  X  t0  ♦  f(n  +  ^ 

Cs 


t 


O^)  For  s  ■  0,1,..., N  ,  add 

r 

r 

♦8<n)  *  -j-  to  $  r  (n  +  1)  . 

Cs 

From  (J»(n)  ,  we  can  obtain  ir(t)  rather  easily  using  equation 
(1),  and  use  ir(t)  to  obtain  system  performance  measures,  such  as  sys- 
tem  availability,  which  we  define  as  the  probability  that  the  desired 
number  of  components  (M)  are  operating  at  time  t  ;  that  is. 


Y 

A(t)  -  l  tt  (t) 
s-0  8 

Various  other  measures  such  as  expected  backorder  level,  expected 
number  of  units  operating,  etc.,  can  also  be  readily  calculated  from 
ir(t)  . 

Computational  efficiency  for  system  performance  measures  can  be 
gained  by  using  $(n)  and  then  converting  to  continuous  time  (see  Gross 
and  Miller,  op.  oit.),  for  example,  defining 
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An  -  I  ♦,<»>  . 
s*0 


we  have 


T  “lit  /  »  \  41 

A(t)  ■  Jo  A»  ■ 


-At. 


This  algorithm  has  been  coded  as  an  interactive  FORTRAN  program 
called  REPTRAN1.  Sample  input  and  output  [A(t)  vs.  t]  are  shown  in  Fig¬ 
ures  2  and  3,  respectively.  Also  shown  is  the  equivalent  (1,1, 1/1,0/%*) 
calculations— the  Dyna-METRIC  model.  This  has  also  been  coded  up  as 
part  of  the  REPTRAN1  program. 

Finite  source,  finite  repair  models  (-,-,-/-,-/f ,f )  are  referred 
to  as  "closed  queuing  network"  models,  for  these  types  of  models  can 
indeed  be  viewed  as  closed  queuing  networks  [the  (1,1, 1/1,0/*,*)  of 
REPTRAN1  is  the  simplest  closed  network,  namely,  a  two-stage  cyclic 
queue] . 


5.  SPECIAL  CASE:  (2,2,2/1,0/f ,f ) 

We  now  present  an  implementation  of  the  techniques  discussed  pre¬ 
viously  to  the  computation  of  transient  probabilities  and  availabilities 
of  a  (2,2,2/1,0/f ,f)  system.  (It  has  been  coded  as  a  FORTRAN  program 
called  REPTRAN2.)  The  system  is  shown  in  Figure  4.  Items  of  one  type 
move  around  a  network  with  six  nodes,  namely,  operational  at  Base  1  (BUI), 
in  or  awaiting  repair  at  Base  1  (BR1) ,  operational  at  Base  2  (BU2) ,  in  or 
awaiting  repair  at  Base  2  (BR2),  operational  (ready  spares)  at  depot  (DU), 
and  in  or  awaiting  repair  at  depot  (DR).  The  description  of  the  system 
will  follow  the  approach  used  earlier.  The  number  of  operating  machines 
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♦  RUN  REPTRAN1 


THIS  IS  REPTRAN 1  PROGRAM 


DO  YOU  WANT  A  HARDCOPY  (Y  OR  N)  ? 

N 

DO  YOU  UANT  TO  PRINT  ALL  OUTPUT  <Y  OR  N)  ? 

Y 

Types  Initial  tine  (assume  zero), final  tine,  tine  increment 
0,15,1 

NUMBER  OF  TIME  POINTS  *  16  NUMBER  OF  TIME  PERIODS  =  15 

Case  nunber:  1 

DO  YOU  UANT  TO  RUN  A  CLOSED  QUEUEING  NETUORK  OR  A  DYNAMETRIC  MODEL 
TYPE  1  OR  2  ACCORDINGLY 
1 

Types  M,C,Y, EPSILON 
3,2,2,-001 

Type  nunber  of  lambda  values  to  be  used 

<7 

4m 

Type  lanbda  values 

*7  “Z 

m  4m  f  m  w 

Type  tines  of  shift  of  lanbdas  (start  w/  0) 

0,6 

Type  nunber  of  nu  values  to  be  used 
2 

Type  nu  values 
.5,-73 

Type  tines  of  shift  of  nu  (start  w/0) 

0,10 

Do  you  want  to  type  in  an  initial  prob.  vector  ?  (Type  Y  or  N) 

If  you  type  N,  program  assumes  ( 1  , 0 , 0 ,  -  -  -  , 0 )  as 
initial  prob.  vector. 

N 


DO  YOU  UANT  TO  RUN  ANOTHER  CASE? 

Y 

Case  numbers  2 

DO  YOU  UANT  TO  RUN  A  CLOSED  OUEUEING  NETUORK  OR  A  DYNAMETRIC  MODEL 
TYPE  1  OR  2  ACCORDINGLY 
2 

TYPE  LAMBDA1 ,LAMBDA2  AND  THE  TIME  AT  UHICH  THE  CHANGE  OF  LAMBDA  OCCURS 

. 6 , .9 ,6 

TYPE  MU  1 , MU2 , MU3 
.5  .5  . 75 

TYPE  THE  TIMES  UHERE  THE  MUS  CHANGE  VALUES 

Figure  2. — Sample  input  for  REPTRAN1:  Closed  queuing  network 
and  Dyna-METRIC  runs. 
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AVAILABILITY 


AVAILABILITY  VS.  TIME 


1 


2 


2 


t 


2 


1 


1 


1 


1 


* 

1 

* 

# 

* 

0.3  +***-»*****  +  *#ii*  +  *********  +  ***4!4-*#**  +  ****  + ****♦**♦*  + ****  +  ****  +  ****  +  :»:  ***  +  *; 

0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0  11.0  12.0  13.0  14.0 

0 


TIME 

Figure  3. — Sample  output  for  REPTRAN1:  Closed  queuing  network 
(1)  and  Dyna-METRIC  (2)  runs. 
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Base  X 


Figure  4. — General  schematic  for  a 
(2,2,2/1,0/f ,f)  repair¬ 
able  item  system. 


at  Base  1  equals  the  maximum  of  the  desired  number,  MSI,  and  the  number  of 
machines  at  node  BUI  (#BU1);  similarly  for  Base  2.*  The  number  of  spares 
available  at  the  depot  spares  pool  is  the  number  at  node  DU  (#DU) .  The  num¬ 
ber  of  busy  repair  channels  at  Base  1  equals  the  minimum  of  the  number  of  re¬ 
pair  channels,  BC1,  and  the  number  of  items  at  node  BR1  (#BR1);  similarly  for 
Base  2  and  the  depot.  Thus,  the  system  parameters  (failure  and  repair  rates) 
and  the  number  of  machines  at  each  node  give  us  total  information  about  the 
system.  Knowledge  of  how  machines  move  around  the  network  completes  the  de¬ 
scription  of  the  system.  The  only  point  which  must  yet  be  specified  in  detail 
is  the  assignment  rule  for  filling  backorders  to  the  depot. 

Consider  the  situation  where  the  number  of  machines  at  node  DR  (#DR) 
is  greater  than  the  spares  allocated  to  the  depot,  DS.  In  this  case  there 
will  be  no  machines  at  node  DU  and  a  backorder  level  of  #DR  -  DS  at  the  depot. 


*Quantities  with  "#"  preceding  them  are  system  state  variables;  quanti¬ 
ties  without  #  are  either  node  designators,  event  descriptors,  or  preset  param¬ 
eter  values. 
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The  question  arises  concerning  the  allocation  of  the  next  aachine  to  complete 
repair  at  the  depot:  To  which  base  should  it  be  given?  We  oust  define  an  al¬ 
location  function  based  on  the  state  of  the  system.  For  now,  let  us  define 
the  state  of  the  system  as  the  number  of  machines  at  each  node*: 

S  -  (#BU1,  #BR1,  #BU2,  #BR2,  #DR,  #DU)  . 

The  number  of  backorders  at  the  depot  from  Bases  1  and  2,  respectively,  are 

#BD1  -  BS1  -  (#BU1  +  #BR1) 

#BD2  -  BS2  -  (#BU2  +  #BR2) 

where  BS1  and  BS2  are  the  numbers  of  spares  allocated  to  Bases  1  and  2, 
respectively.  Possible  allocation  functions  are 

{1  if  #BD1  >  #BD2 

2  if  #BD2  >  #BD1 

or 

{1  if  #BD1/BS1  >  # BD2/BS2 

2  if  #BD1/BS1  >  0BD2/BS2 

In  the  program  REPTRAN2  we  use  a  randomized  generalization  of  the  second  alloca¬ 
tion  function:  Let  PROB.ALL(*>  be  a  function  that  equals  the  probability  of 
assigning  the  repaired  machine  to  Base  1.  The  assignment  rule  is  given  in 
Table  II.  The  weights  W1  and  W2  are  supplied  by  the  user  and  reflect  his 

TABLE  II 

ASSIGNMENT  RULE  FOR  SENDING  A  REPAIRED  MACHINE 
FROM  DEPOT  IN  BACKORDER  SITUATION 


Condition 

Probability 
of  Assignment 
to  Base  1 

Probability 
of  Assignment 
to  Base  2 

W1  *  #BD1  >  W2  * 

#BD2 

1 

0 

W1  *  #BD1  -  W2  * 

#BD2 

1/2 

1/2 

W1  *  #BD1  <  W2  * 

#BD2 

0 

1 

*Because  it  is  a  closed  system,  we  shall  see  later  that  £  can  be 
characterized  by  fewer  state  variables. 
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strategy  for  favoring  one  base  over  the  other  when  both  have  backordered 
units. 

The  description  of  the  system  is  now  complete  and  we  proceed  to 
describe  the  Markov  model  of  the  system  using  the  SERT  approach.  This 
requires  (i)  identifying  the  state  space  (as  a  vector),  (li)  defining 
the  event  set,  and  for  each  event  (ill)  computing  the  vector  of  transi¬ 
tion  rates,  and  (iv)  computing  the  vector  of  target  states.  The  randomi¬ 
zation  algorithm  can  then  be  used  to  compute  transient  probabilities  to 
any  user-specified  accuracy  (e) .  The  events  are  given  in  Table  III; 
once  the  state  space  is  described  as  a  one-dimensional  vector,  steps 
(iii)  and  (iv)  are  straightforward. 

In  general,  the  state  space  appears  to  have  six  dimensions,  but 
because  of  one-for-one  ordering  and  conservation  of  the  total  number  of 
items  in  the  system,  the  state  space  actually  has  a  lower  dimension. 

TABLE  III 

THE  EVENTS  THAT  ACCOUNT  FOR  ALL  THE  STATE  CHANGES  OF 
A  (2,2,2/1,0/f ,f)  SYSTEM  WITH  A  SINGLE  TYPE  OF  ITEM 

Name  Description 

B1R  Repair  completed  at  Base  1 

B2R  Repair  completed  at  Base  2 

FIB  Failure  at  Base  1  (base  repairable) 

F2B  Failure  at  Base  2  (base  repairable) 

FID  Failure  at  Base  1  (depot  repairable) 

F2D  Failure  at  Base  2  (depot  repairable) 

DR1  Repair  completed  at  depot  and  sent  to  depot 
spares  pool  if  not  backorder  situation; 
otherwise  sent  to  Base  1 

DR2  Repair  completed  at  depot  while  backorder 
situation  and  sent  to  Base  2 
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The  description  of  the  state  space  breaks  into  two  situations:  no  depot 
spares  available,  and  some  depot  spares  available, 

SmSoKls* 

where 

Sq  -  states  with  depleted  depot  spares  pool 
S+  *  states  with  nondepleted  depot  spares  pool. 

First  consider  Sg  .  In  this  case,  it  is  possible  to  describe  the  state 
of  the  system  with  four  numbers: 

(#BU1,  #BD1,  //BU2,  #BD2) 

It  is  known  that  it  DU  *  0  and  the  remaining  machines  are  at  node  DR. 

The  feasible  states  of  Sq  are  subject  to  two  constraints: 

#BU1  +  #BD1  <  BS1 
#BU2  +  it BD2  <  BS2 


and  thus  SQ  is  a  Cartesian  product. 


VT2 


where 


T  -  {(0BD1,  0BU1) :  //BD1  +  //BUI  <  BSl) 

T2  -  {<0BD2,  f/BU2) :  //BD2  +  0BU2  <  BS2) 

These  sets  are  shown  in  Figure  5.  (The  notation  T  is  used  because  the 
spaces  are  triangular.)  Note  that  the  number  of  points  in  and  T ^  are 

IT  I  .  (BS1+1) (BS1+2) 

1  1'  2 

and 

|T  |  .  (BS2+mBS2+2) 

respectively,  and  the  number  of  states  in  SQ  is  the  product 

(5  |  »  (BS1+1)  (BS1+2)  (BS2+1)  (BS2+2) 
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#BD1  2  • 


#BU2 

0  12  3 


4  5 


Figure  5. — Examples  of  state  space 
to  describe  individual 
bases;  in  the  case  BS1  - 
4  and  BS2  -  5. 


Now  let  us  consider  the  states  where  the  spares  pool  at  the  depot  is  not 
empty,  S+  .  In  this  case  the  state  of  the  system  can  be  described  by 


three  numbers: 


<#BU1,  #BU2,  #DU) 


The  constraints  on  these  are 


#BU1  <  BS1 
#BU2  <  BS2 


1  <  #DU  <  DS 


We  can  condition  on  the  value  of  #DU  to  get  S .  into  the  form 
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S  KJ  S  \J  S  Kl  S 

>  V  V  •  “  *DS 

where  consists  of  states  with  exactly  i  machines  in  the  depot 
spares  pool;  such  a  set  is  depicted  in  Figure  6.  Note  that  each  is  a 


rectangle  and 


|S  |  -  (BS1  +  1)  (BS2  +  1) 


(BS1  +  1) (BS2  +  1)DS 


and  the  total  number  of  states  is 


isi  -  +  <BS1  +  1KBS2  +  i)DS 


Examples  of  state  space  sizes  are  given  in  Table  IV. 

Although  in  principle  we  can  apply  the  randomization  algorithm  to 
a  system  whose  state  space  is  described  in  this  complex  multidimensional 
form,  we  choose  to  work  with  one-dimensional  state  spaces.  We  made  this 
choice  for  two  reasons.  First,  one  general  implementation  of  the  ran¬ 
domization  algorithm  will  work  on  all  systems  after  they  are  put  into 


#UB2 

0  1  2  3  4  5 

o  •  •  •  •  •  • 


^  •  •  •  •  •  • 

#UB1  2 . S1 

3  *****  # 

4  •  •  •  •  •  • 


Figure  6. — The  state  space  describing 
status  of  bases  given  that 
depot  spares  are  available; 
BS1  -  4  and  BS2  -  5. 
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TABLE  IV 

SIZE  OF  THE  STATE  SPACES  OF  (2,2,2/1,0/f ,f)  SYSTEM 
FOR  SELECTED  VALUES  OF  BS1,  BS2,  AND  DS 


BS1 

BS2 

DS 

l$0l  + 

isj 

«  |S| 

2 

2 

2 

36 

18 

54 

4 

4 

2 

225 

50 

275 

6 

6 

2 

784 

98 

882 

8 

8 

2 

2025 

162 

2187 

10 

10 

2 

4356 

242 

4598 

12 

12 

2 

8281 

338 

8619 

18 

18 

2 

36100 

722 

36822 

24 

24 

2 

105625 

1250 

106875 

one-dimensional  form,  and  it  seems  to  decrease  the  problems  in  verifying 
the  program.  Second,  and  even  more  important,  working  directly  with 
vectors  rather  than  structures  such  as 


DS 

S  m  T,  x  TjV  Si 
i-1 

will  speed  up  the  algorithm  by  a  significant  factor  because  it  is  much 
faster  to  "fetch"  and  "store"  an  element  of  a  one-dimensional  vector 
than  it  is  to  "fetch"  and  "store"  an  element  of  a  complicated  multidi¬ 
mensional  array.  We  estimate  that  this  speeds  up  REPTRAN2  by  a  factor 
of  four  or  more. 

Thus,  let  us  consider  putting  the  elements  of 

DS 

5  -  T  x  T 2KJ  \J  S 
L  i-1  1 


-  20  - 


T-468 


tnco  a  one-dimensional  form.  First  consider  T^;  this  triangular  shaped 
region  can  be  transformed  into  a  linear  space  i.^  by  placing  the  rows  end 
to  end, 

Lx  -  (00,  01,  02,  03,  04,  10,  11,  12,  13,  20,  21,  22,  30,  31,  40} 
Doing  the  same  thing  for  T ^  yields  Lj »  thus, 

S0  *  T1XT2  2  L1XL2 

becomes  a  two-dimensional  set  which  can  be  linearized  again  by  placing 
the  rows  end  to  end  as  shown  graphically  in  Figure  7.  The  first  |Sq| 
elements  of  our  one-dimensional  vector  will  correspond  to  the  states  in 
Sq  in  linear  order  as  presented  in  Table  IV. 

The  remaining  |S+|  elements  in  the  one-dimensional  vector  will 
correspond  to  the  states  in  S+  arranged  linearly  as  follows: 

DS 

-  y  si  • 

i*l 

where  is  depicted  in  Figure  6.  Each  is  put  into  one-dimensional 
form  by  putting  rows  end  to  end  (as  in  Figure  7).  Finally,  these  DS 


50  "  L1  X  L2 


1  02 


2  dimensional 
version 


00  01  02  50 

L2 

0000,  0001 .  0050,  0100,  0101,  ....  0150,  ....  4000,  4001,  ....  4050 

Linear  version 

Figure  7. — Depiction  of  linearization  of  the  Cartesian  product 
being  transformed  into  a  one-dimensional  space. 


JK  iVjf!'-*-  :<■'  . 


-  21  - 


T-468 


one-dimensional  S^'s  are  put  end  to  end  starting  with  and  ending  with 
(Thus  the  last  state  in  this  one-dimensional  listing  is  (BS1,  BS2, 

DS) ,  which  corresponds  to  BS1  units  "up"  at  Base  1  (#BU1  *  BS1) ,  BS2  units 
"up"  at  Base  2,  and  OS  units  in  the  spares  pool  at  the  depot,  i.e.,  a 
perfect  system  with  no  failed  units.  The  first  state  listed  in  the  one- 
dimensional  version  of  the  state  space  is  (0,  0,  0,  0),  which  corresponds 
to  BS1  units  in  repair  at  Base  1,  BS2  units  in  repair  at  Base  2,  and  DS 
units  in  repair  at  the  depot.) 

This  completes  the  description  of  the  state  space  of  this  (2,2,2/ 
l,0/f,f)  system  as  a  one-dimensional  vector.  It  is  used  by  the  program 
REPTRAN2 .  A  major  part  of  REPTRAN2  is  subroutines  that  take  the  param¬ 
eters  and  rates  and  compute  one-dimensional  vectors  of  transition  rates 
and  one-dimensional  vectors  of  target  states  for  each  event  in  Table  III. 

It  then  uses  the  randomization  algorithm  logic  to  compute  transient  state 
probabilities  and  availabilities  as  described  in  Section  2. 

The  program  REPTRAN2  can  also  compute  transient  probabilities  for 
(2, 2,2/1,0/f ,f )  systems  whose  underlying  rates  (failure  and  repair  rates) 
change  at  discrete  points  in  time  (on  a  lattice  time  scale) — up  to  five 
changes  are  accommodated  for  each  underlying  rate. 

We  note  that  by  setting  certain  parameters  to  special  values, 

REPTRAN2  can  handle  the  corresponding  special  cases.  If  DS  *  0,  we  lose 
the  second  level  of  supply  and  get  a  (2,2,1/1,0/f ,f)  system.  If  FBI  ■  1 
then  all  failures  at  Base  1  are  base  repairable  and  Base  1  does  not  in¬ 
teract  with  the  depot:  thus  Base  1  is  a  (1,1,1/1,0/f ,f)  system  and  Base 
2  and  the  depot  form  a  (1, 2,2/1,0/f ,f)  system,  or  (1,2,1/1,0/f ,f )  if  DS  ■  0. 


-  22  - 


T-468 


6.  COMPUTER  RUNS 

The  REPTRAN2  program  has  been  run  for  various  (2,2,2/1,0/f ,£) 
systems  as  described  in  Section  5.  The  program  executes  In  quite  short 
times  for  systems  with  large  state  spaces.  It  also  compares  the  prob¬ 
abilities  computed  by  the  Dyna-METRIC  model  for  corresponding 
-,-/00,00)  systems,  which  are  also  a  part  of  the  REPTRAN2  code.  [For  the 
Dyna-METRIC  developments,  details  for  these  corresponding  systems,  see 
Gross,  Kioussis,  and  Miller  (1982).]  Furthermore,  it  is  reassuring  that 
when  REPTRAN2  was  run  on  a  system  whose  parameters  gave  rise  to  a 
(1,1,1/1,0/f ,f)  subsystem  (as  discussed  in  Section  4),  the  program  gave 
the  same  answers  as  REPTRAN1. 

REPTRAN2  is  interactive  and  the  user  has  the  option  of  running 
either  the  Dyna-METRIC  model  or  the  closed  queuing  network  model  for  a 
(2, 2, 2/1,0/* , •)  system.  Availability  is  defined  as  the  probability  that 
the  number  of  machines  "up"  at  a  base  meets  or  exceeds  the  number  desired 
to  be  operating,  and  is  the  primary  measure  outputted. 

We  ran  a  few  sample  cases  to  get  an  idea  of  how  well  this  program 
performs  (with  reference  to  time)  in  computing  exact  transient  solutions. 
The  sample  cases  we  ran  are  described  in  Table  V.  For  each  case,  all  the 
necessary  parameters  (which  the  user  must  supply)  are  listed.  The  fail¬ 
ure  rates  and  repair  rates  shift  at  the  times  indicated  for  each  case. 


For  example,  in  Case  la:  At  time  tQ  -  0  ,  the  failure  rate  of  a 
machine  at  Base  1  is  A.  ■  .4  ;  then  at  time  t,  ■  6  ,  this  rate  changes 
to  A^  ■  .6  .  At  time  tg  ■  0  the  repair  rate  of  a  single  repair  chan¬ 


nel  at  Base  1  is 


then  at  time 


“  10  ,  this  rate  changes 


to  *  .75  .  Similar  changes  occur  for  the  failure  and  repair  rates 


at  Base  2  and  the  repair  rate  at  the  depot. 
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TABLE  V 

PARAMETER  VALUES  OF  THE  SAMPLE  CASES  SOLVED  USING  REPTRAN2 

A.  Input  Parameter  Names 
Symbol  Meaning 


BSi 

MSi 

BCi 

FBI 

C0 

tXi 

w 

Pi(t0) 

Ml<tU1) 

Wi 

DS 

DC 

uD<t0) 


Allocation  of  total  stock  to  Base  i  (operating  machines  plus 
spares),  i  *  1,2 

Desired  number  of  working  machines  at  Base  i 
Number  of  repair  channels  in  repair  shop  at  Base  i 
Proportion  of  items  failing  at  Base  i,  base  repairable 
Time  zero 

Time  of  shift  in  mean  failure  rate.  Base  i 
Time  of  shift  in  mean  repair  rate.  Base  i 
Initial  mean  failure  rate.  Base  i  items 
Shifted  mean  failure  rate.  Base  i  items 
Initial  mean  repair  rate.  Base  i  items 
Shifted  mean  repair  rate.  Base  i  items 

Weighting  factor  for  filling  depot  backorders  for  Base  i 

Number  of  depot  spares 

Number  of  depot  repair  channels 

Initial  mean  depot  repair  rate 

Time  of  shift  in  mean  depot  repair  rate 

Shifted  mean  depot  repair  rate 

Error  tolerance 
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Table  V — continued 


B.  REPTRAN2  Run  Description 

BS1 

MSI 

BC1 

FBI 

w 

xi<V 

ui(V 

BS2 

MS2 

BC2 

FB2 

X2(tX2> 

^2^C0^ 

U2(tU2) 

DS 

DC 

W1 

W2 

w 

^D^D^ 

e 

Cases 

la 

4 

2 

2 

.7 

.4(0) 

.6(6) 

.5(0) 

.75(10) 

5 

3 

2 

.5 

.4(0) 

.6(8) 

.6(0) 

.9(12) 

2 

2 

.4 

.6 

.3(0) 

.45(11) 

.001 

lb 

4 

2 

2 

.7 

.4(0) 

.6(6) 

.5(0) 

.75(10) 

5 

3 

2 

.5 

.4(0) 

.6(8) 

.6(0) 

.9(12) 

2 

2 

.5 

.5 

.3(0) 

.45(11) 

.001 

2a 

4 

2 

2 

.7 

.2(0) 

.4(6) 

.5(0) 

.75(10) 

5 

3 

2 

.5 

.1(0) 

.15(8) 

.4(0) 

.6(12) 

2 

2 

.4 

.6 

.3(0) 

.45(11) 

.001 
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The  seven  cases  in  Table  V  were  solved  exactly  over  16  time 
points  {0,1, . . . ,15}.  The  CPU  execution  times  are  shown  in  Table  VI,  as 
well  as  the  size  of  the  state  space  of  the  model.  Also  shown  is  the  set 
of  time  points  at  which  rate  changes  (failure  or  repair)  occurred. 
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TABLE  VI 

SIZE  OF  STATE  SPACE  AND  RUNNING  TIMES  FOR  REPTRAN2 
SOLUTION  OF  VARIOUS  (2,2,2/1,0/f ,f)  SYSTEMS 


System 

Number 

Size  of 
State  Space 

Homogeneous  Time 
Interval 

CPU 

(sec.)a 

la 

375 

0,  6, 

8, 

10, 

11. 

12 

10.18 

lb 

375 

0,  6, 

8, 

10, 

11. 

10 

10.48 

2a 

375 

0,  6, 

8, 

10, 

11. 

12 

8.07 

2b 

375 

0,  6, 

8, 

10, 

11. 

12 

8.87 

3 

375 

0,  6, 

8, 

10, 

11. 

12 

9.04 

4 

3366 

0,  6, 

8, 

10, 

11, 

12 

190.25 

5 

20748 

0,  6, 

10 

1393.73 

aRunning  time  is  approximately  proportional  to  size  of 
state  space  x  number  of  event  types  x  number  of  occur- 
rences  of  Poisson  process. 


The  seven  cases  in  Table  VI  were  also  solved  using  the  Dyna- 
METRIC  approximate  model  (2,2,2/l,0/»,») .  For  this  comparison  it  was 
necessary  to  select  failure  arrival  rates  at  each  base.  We  chose  an  ar¬ 
rival  rate  of  MSI  *  at  the  Base  1  repair  shop  and  MS2  *  at  the 
Base  2  repair  shop.  For  systems  operating  at  moderate  to  high  availa¬ 
bilities,  these  should  be  approximately  correct. 

We  show  the  computer  input  and  output  for  case  5  given  in  Table  V 
as  Figures  8  through  12.  The  Interactive  input  is  shown,  first  modelled 
as  an  exact  queuing  network  (2,2,2/1,0/f ,f)  and  then  modelled  as  a  Dyna- 
METRIC  model  (2, 2, 2/1,0 /•,»).  Both  models  were  solved  for  availabilities 
and  then  superimposed  on  three  plots.  The  exact  availabilities  are 
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TYPE:  INITIAL  TIME  (assune  zero), FINAL  TIME,  TIME  INCREMENT 
0,15,1 

NUMBER  OF  TIME  POINTS  =  16  NUMBER  OF  TIME  PERIODS  *  15 

CASE  NUMBER  :  1 

DO  YOU  UANT  TO  RUN  A  CLOSED  QUEUEING  NETWORK  OR  A  DYNAMETRIC  MODEL 
TYPE  1  OR  2  ACCORDINGLY 
1 

TYPE  BS1 ,MS1 ,BC1 , AND  FBI 

18.14.2, -6667 

BS 1  :  IS  MSI  :  14  BC1  :  2  FBI  sO.6667 

TYPE  BS2,MS2,BC2,AND  FB2 

13.10.2, -6667 

BS2  :  13  MS2  :  10  '  BC2  s  2  FB2  .*0.6667 

TYPE  DS,0C,U1 , U2 
3,4, .5, .5 

DS  :  3  DC  :  4  U1  :0.5000  142  :0.5000 

TYPE  NUMBER  OF  FAILURE  RATES  TO  BE  USED  IN  BASE  1 

2 

TYPE  THE  FAILURE  RATES 
-2,  .3 

TYPE  THE  TIME  UHERE  THOSE  CHANGES  IN  RATES  OCCUR  (START  UITH  0) 

0 , 6 

TYPE  NUMBER  OF  REPAIR  RATES  TO  BE  USED  IN  BASE  1 

n 

4b 

„.JYP.c  the  .REPAIR  RATES 

1  -  ,  1  .5 

TYPE  THE  TIME  UHERE  THOSE  CHANGES  IN  RATES  OCCUR  (START  UITH  0) 

0.10 

TYPE  NUMBER  OF  FAILURE  RATES  TO  BE  USED  IN  BASE  2 

2 

TYPE  THE  FAILURE  RATES <  BASE  2) 

-143,-2143 

TYPE  THE  TIME  UHERE  THOSE  CHANGES  IN  RATES  OCCUR  (START  UITH  0) 

0,6 

TYPE  NUMBER  OF  REPAIR  RATES  TO  BE  USED  IN  BASE  2 

2 

TYPE  THE  REPAIR  RATES  (BASE  2) 

1  -0, 1  .5 

TYPE  THE  TIME  UHERE  THOSE  CHANGES  IN  RATES  OCCUR  (START  WITH  0) 

0,10 

TYPE  NUMBER  OF  REPAIR  RATES  TO  BE  USED  IN  DEPOT 

2 

TYPE  THE  REPAIR  RATES  (DEPOT) 

-5,-75 

TYPE  THE  TIME  UHERE  THOSE  CHANGES  IN  RATES  OCCUR  (START  UITH  0) 

0,10 

TYPE  THE  DESIRED  MAGNITUDE  OF  THE  ERROR ( EPS  I  LON ) 

.001 

Figure  8. — Sample  input  for  REPTRAN2  case  5,  closed  queuing  network. 
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00  YOU  UANT  TO  RUN  ANOTHER  CASE?(Y  or  N> 

Y 

CASE  NUMBER  :  2 

DO  YOU  UANT  TO  RUN  A  CLOSED  QUEUEING  NETUORK  OR  A  DYNAMETRIC  MODEL 
TYPE  1  OR  2  ACCORDINGLY 
2 

TYPE  THE  TIME  AT  UHICH  THE  FAILURE  RATE  AND  THE  REPAIR  RATE  CHANGE  AT  BAS! 

6,10 

TYPE  THE  TUO  FAILURE  RATES  FIRST  AND  THEN  THE  TUO  REPAIR  RATES 
2.0  4.2  1.0  1.5 

TYPETHESTOCK  LEVEL  AT  BASE  1  AND  THE  PROPORTION  OF  THE  FAILED  ITEMS  G0I1 
BASE  REPAIR  1 

4. . 6667 

TYPE  THE  TIME  AT  UHICH  THE  FAILURE  RATE  AND  THE  REPAIR  RATE  CHANGE  AT  BAS! 

6,10 

TYPE  THE  TUO  FAILURE  RATES  FIRST  AND  THEN  THE  TUO  REPAIR  RATES  (BASE  2) 
1.43,2.143,1.0,1.5 

TYPE  THE  STOCK  LEVEL  AT  BASE  2  AND  THE  PROPORTION  OF  THE  FAILED  ITEMS  GOIf 
CASE  REPAIR  2 

3. . 6667 

TYPE  THE  TIME  AT  UHICH  THE  REPAIR  RATE  CHANGES  AT  THE  DEPOT 
10 

TYPE  THE  TUO  REPAIR  RATES  (FOR  DEPOT) 

.5, .75 

TYPE  THE  STOCK  LEVEL  AT  THE  DEPOT 
3 

Figure  9. — Sample  input  for  REPTRAN2  case  5,  Dyna-METRIC .  J 


plotted  using  the  symbol  "1"  and  the  approximate  Dyna-METRIC  availabili¬ 
ties  are  plotted  using  the  symbol  "2."  For  time  points  where  the  symbol 
"1"  fails  to  appear,  it  coincides  with  the  "2."  "Availability  1"  is 
availability  at  Base  1,  "Availability  2"  is  availability  at  Base  2,  and 
"Availability  3"  is  simultaneous  availability  at  both  Base  1  and  Base  2. 
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AVAILABILITY  t 

1.0  2  2 

*  1  2 

*  .  ....  .  .  _  2 

*  \ 

* 

* 

* 

*  1 

0.8  ♦ 

* 

* 

* 

* 

* 

* 

* 

0.6  + 

* 

* 

* 

* 

* 

* 

* 

0.4  + 

* 

* 

* 

# 

+ 

* 

+ 

0.2  + 

* 

+• 

* 


AVAILABILITY  1  VS.  TIHE 


2 


2 


2 


2 


2 


1 


1 


1 


2 


1 


1 


1 


1 


1 


* 


0.0  +  * * * 

0.0  1.0 

0 


2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0  11.0  12.0  13.0  1 

TIME 


Figure  10. — Sample  output  for  REPTRAN2  case  5, 
1:  closed  queuing  network, 

2:  Dyna-METRIC, 
availabilities  at  Base  1. 
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AVAILABILITY 
2 
1 

•8 

* 

* 

* 

O.Q  > 


2  2 


AVAILABILITY  3  VS.  TIME 


2  2 


v .  •> 


♦ 

* 

* 

Us 

* 

* 

+ 

4s 

+ 

* 

* 

* 

* 


1  2 
1  2 


0, 


* 

* 

* 

n 

* 

* 

j. 

* 


:8 
* 

0.0  +  sfe  A  ♦  s8  +♦♦♦♦  +  ♦**♦  +  ♦♦♦♦♦♦♦♦♦  +  *♦♦♦♦♦♦♦♦♦*♦♦#  +  ♦♦♦#♦♦♦♦♦♦♦**  s|i  ♦♦♦*♦♦♦♦♦#♦♦  •■ 
0.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0  11.0  12.0  13.0 

0 

TIME 

Figure  12.— Sample  output  for  REPTBAN2  case  5, 

1:  closed  queuing  network, 

2:  Dyna-METRIC , 

simultaneous  availabilities  at 
both  bases. 
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7 .  CONCLUSIONS 

The  SECT  modelling  technique  and  randomization  computing  algo¬ 
rithm  [methodology  developed  in  Gross  and  Miller  (1982)]  has  been  ap¬ 
plied  and  implemented  for  computing  transient  performance  measures  of 
multi-echelon  repairable  item  inventory  systems.  We  have  shown  that  it 
is  feasible  to  compute  exact  probabilities  for  systems  with  large  state 
spaces  (20,000  or  more  states).  Furthermore,  for  the  types  of  systems 
under  consideration,  we  believe  that  significantly  larger  cases  are 
feasible  using  a  truncated  state  space  approach;  that  is,  lumping  the 
vast  number  of  very  low  probability  states  together  as  one. 

For  example,  in  a  (1,1,1/1,0/f ,f)  case  with  31  states  (M  »  20, 

Y  ■  10),  we  found  that  the  probabilities  of  s  units  in  resupply  for 
s  >  15  was  zero  to  at  least  three  significance  figures.  Lumping  states 
s  «  15,  16,  ...,  31  together  would  reduce  the  problem  from  31  states  to 
16,  a  savings  of  almost  50%.  We  estimate  in  the  20,000  state  space 
example  that  such  a  procedure  would  easily  cut  the  number  of  states  in 
half. 

We  have  used  the  SECT  modelling  technique  on  a  (b,2,2/k,l/f ,f ) 
system  [see  Gross,  Kioussis,  and  Miller  (1982)].  This  has  not  been 
coded,  but  for  moderate  b  and  k  ,  using  the  truncated  state  space 
approach,  development  of  an  efficient  code  should  be  feasible.  Con¬ 
ceptually,  of  course,  the  most  general  (b,r,s/k,j/f ,f )  system  could  be 
modelled  using  SECT;  the  problem,  of  course,  is  the  state  space  size  for 
cases  other  than  those  with  very  small  values  of  b  ,  r  ,  s  ,  k  ,  and  j 
The  truncated  state  space  approach  offers  the  most  promise  for  treating 


these  models. 
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